Detecting the concurrence of an unknown state with a single observable 
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While the detection of entanglement has been proved already to be quite a difficult task, ex- 
perimental quantification of entanglement is even more challenging. In this work, we derive an 
analytical lower bound for the concurrence of a bipartite mix quantum state in arbitrary dimension. 
The lower bound is experimentally implementable in a feasible way, which enabling quantification 
of entanglement in a broad variety of cases. 
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Entanglement is a distinctive feature of quantum me- 
chanics 0, and an indispensable ingredient in vari- 
ous kinds of quantum information processing applications 
vary from quantum cryptography 0] and quantum tele- 
p_ortation |J| to measurement-based quantum computing 

The use of entanglement as a resource not only bears 
the question of how it can be detected, but also how it 
can be quantified. For this purpose, several entangle- 
ment measures have been introduced, one of the most 
prominent of which is the concurrence Q . 

However, calculation of the concurrence is a formidable 
task as the Hilbert space dimension is increasing. Good 
algorithms and progresses have been obtained concern- 
ing lower bounds [9l4l3|. Considerable progress is made 
in order to give a pure ly alge braic lower bound and ex- 
perimental verifying it 131-4231] . 

We start with a generalized definition of concurrence 
for a pure state \ijj) in the tensor product Ha'&'Hb of two 
Hilbert spaces Ha,T~Ib for systems A, B.The concurrence 
is defined by 
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where the reduced density matrix pa is obtained by trac- 
ing over the subsystem -B.The concurrence is then ex- 
tended to mixed states p by the convex roof, 
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for all possible ensemble realizations p = 'YlnPi^i) 
where pi > and '^2 i Pi — 1. From definition, a state p 
is separable if and only if C(p) = 0. 

In this letter, we will show that strong lower bounds of 
concurrence can be derived by exploiting close relations 
between concurrence and a recently introduced detection 
criteria for bipartite entanglement. The lower bound we 
found is not only analytical, but also is physical accessi- 
ble, i.e., we can easily use the exis ting exper iments tech- 
nology to obtain this bound(e.g., |l7l4l9l l25l]). 

To disclose the connection between concurrence and 
entanglement detection criteria, let us first review the 



known critcrions. Until now, many possible ways to 
detect entanglement have been proposed. These range 
from Bell inequalities, entanglement witnesses and spin 
squeezing inequalities to entropic inequalities, the mea- 
surement of nonlinear properties of the quantum state 
and the approximation of positive maps, see[l|, Q. 

Among so many methods, the PPT criterion and the 
realignment criterion are the two distinguished ones, they 
are most powerful and widely used by quantum informa- 
tion community. The Peres-Horodecki criterion of posi- 
tivity under partial transpose (PPT criterion) [2(| [27| say 
that p TA > should be satisfied for a separable state, 
where p TA stands for a partial transpose with respect to 
the subsystem A. This criterion is even to be sufficient 
for 2x2 and 2x3 bipartite systems. However, for high 
dimension system, there exists entanglement states that 
are positive under partial transpose, the so called bound 
entanglement states (BES) [11]. 

Another complementary operational criterion for sep- 
arability called the realignment criterion is very strong 
in detecting many of BES [29| . This criterion states that 
a realigned version IZ(p) of p should satisfy ||7£(p)|| < 1 
for any separable state p. Note that realignment crite- 
rion is in some sense " dual " to the PPT criterion, in 
the following way: A density operator can be written as 



p = J^{ij\p\kl)\ij){M\ 
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A partial transpose with respect to the first system A is 



Pkj,il- 



While the realignment of p is defined as: 
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Recently, a criterion was introduced, which was strictly 
stronger than the realignment criterion. It says that, if a 
state is separable, then the following must holdplol]: 



\\K(p -pa® PB )\\ < ^/(1-Trpi)(l-Trp|)] (6) 
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We now derive the main result of this Letter. 
Theorem: For any m ® n (m < n) mixed quantum 
state p, the concurrence C(p) satisfies the following: 
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Proof. — Without loss of generality, we suppose that a 
pure m ® n (m < n) quantum state p := IV'XV'I has the 
standard Schmidt form 



(8) 



where ^fpl [i = l,...m) are the Schmidt coefficients, 
\di) and \bi) arc orthonormal basis in Ha and Hb, re- 
spectively. The two reduced density matrices pa and pb 
have the same eigenvalues of pi. It follows 
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which varies smoothly from 0, for pure product states, to 
2(m — for maximally entangled pure states. 

Then we get 

\\R{p- p A ® pb)\\ 
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So from Cauchy-Schwarz inequality, we get 
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Now, we will show that the inequality (7) also holds 
for mix state. Assume we have found the optimal de- 
composition EiPiPi f° r M to achieve the infimum of 
C(jo), where pi are pure state density matrices. Then 
C(p) = ^2iPiC(pi) by definition. Now, we need to prove 



that after mixture, the bound will not become bigger, 
le.J(J2iPipi) < J2iPif(Pi)- From the convex property 
of the trace norm, we know that ||72-(p — Pa® Mb) II is 
decreasing after convex combination, so the only thing 
left to prove is that y(l — Trp A )(l — Trp 2 B ) is increasing 
after convex combination. It is sufficient to consider the 
case of p = \p\ + \pi- Now, we will prove that 



[l-±Tr( PlA + P2A )2][l- \Tr(p 1B 
>i[l-Tr(^)] + |[l-Tr(^)] 
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Tr(pl B ) := Xl , Tr(p 2 A ) 



For clear, denote Tr(p 2 A ) 
Tr (P2s) : = x 2) Ty(piaP2a) ■= x 3i Tt(pibP2b) ■= %4, then 
the inequality (10) reduced to prove that the function 
F := -^{xi+x 2 ) 2 + \{x 1 +X2){x 3 +x i ) + \x 3 x 4 -\{x 3 + 
X4) + h{%i + X2) > OTo get the minimal value of F, we 
will use the Lagrange multipliers method, and find that F 
attained its minimal value at the point (x\, X2, X3, X4) = 
(1, 1, 1, 1). And in this case, the minimal value is F = 0. 
So we get F > 0. Theorem is proved. ■ 

The most prominent feature of this theorem is that it 
not only allows to obtain a strong lower bound for the 
concurrence without any numerical optimization proce- 
dure, but also this bound is measurable, i.e., it is directly 
accessible in currently existing laboratory experiments. 

We will explain this in detail. We will see that, our 
bound is experimentally implementally by means of local 



observables, use the method as that of [13l.ll8l.l24l.l25l]. 

First, note that the concurrence of a bipartite pure 
state has another representation as 13l. Il8| 



= y/2(l- Trp\) = y/(^\ ® WW) ® W>, (H) 



where A = AP^ ® P_i 2) . 



P_ is the projector on the 



antisymmetric subspace Hi A Hi of the two copies of the 
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Ki = 4(1^) ™ p(2) 
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PT ; ), then 1 - Trp A = ±Tr(p ® pKj), 
^Tr(p ® pK 2 ), so we can obtain the 
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term - Trp 2 A )(l - Trp^l provide that we have two 
copies of the state, see [24|.T25|. 

On the other hand, the term \\TZ(p — pa ® Pb)\\ 
can be obtained by generalized entanglement witness 
27 . 3~0] .Entanglement witnesses (EW) are Hermitian op- 
erators that have positive averages on all separable states, 
but a negative one on at least one entangled state. It was 
shown that, a state is entangled if and only if it is de- 
tected by some EW 27j • EW can be measured locally, 
and one can optimize such measurements in various as- 
pects [32j. Nowadays, entanglement witnesses are rou- 
tinely used in experiments to detect en tang lement (see 
e.g., [33| ). Using the similar method of (34 |. we can di- 
rectly measure the value of \\TZ(p — pa ® Pb)\\ by a gen- 
eralized witness. The method is as follows: from the 
result of matrix analysis [35]. every operator has a singu- 
lar value decomposition (SVD), so we can get the SVD 
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oiRip- p A ®p B ) as R{p~pa®Pb) = UDV + , with U, V 
unitary matrices, D is a diagonal matrix. Then since the 
trace norm has a variational representation as \\A\\ = 
m&x{\Tr[X + AY] |}, where X, Y are unitary matrices, we 
can get that \\K{p- p A ® p B )\\ = Ti-[VU+TZ(p- p A ®p B )}- 
Now the only problem leaved is how to deal with the 
realignment operation. Note that every map <p on m ~ 
ner product space can induce its adjoint map </)*,i.e., 
Tr[<t>{X)Y+] = Tt[X(4>*(Y)) + ], and for our question, it 
reads that Tr[VU + TZ(p - p A ® Pb)} = Tr[K*(VU+)(p - 
Pa® Pb)], and the adjoint map of the realignment oper- 
ation is defined by TZ*(VU+) := [R.- 1 (VU+) T ] T , where 
T is the transpose, and 7?. -1 is the inverse map of the re- 
alignment operation. Now define W := [TZ" 1 (VU + ) T ] T , 
then we get that \\K(p- p A ® p B )\\ = Tr[W(p-_p A ® p B )}, 
which is clearly physical accessible(e.g., see |33(). 

Next we consider some examples to illustrate further 
the tightness and significance of our bound. To show that 
our bound is close to the real concurrence, note that in 
[38| , the authors find a fast optimal algorithm to calculate 
the entanglement of formation of a mixed state, so we can 
use the optimal algorithm of [HI to give the estimation of 
the concurrence, and comparing it with our bound. The 
concurrence using optimal method of [38| is represented 
by blue colors, our low bound is represented by green 
colors. 

Example 1 Isotropic states are a class of U ® U* in- 
variant mixed states in d x d systems 



trix p is real and symmetric, 
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where |*+) = v/V^EiLi \ U ) and F = (*" 
satisfying < F < 1, is the fidelity of pf and 
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(13) 
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where < a < 1. Let us consider a mixture of this state 
with white noise, 



Pip) =PP+0--p)g, 



(14) 



and show the curves 1 
Trp|)/2 - 0, v/(l-Trp;U(l-Tr^) 



r\\ - iTrp 2 A 
\K{p - PA® 



\np)\\ = o, i- 

o%)/2 = 0, V(l - Trpi, 

= with respect to the CCNR criterion, its op- 
timal nonlinear witness, and Theorem 1 in Fig. 2. It is 
found that the state pip) still has entanglement when 
p = 0.9955, a = 0.236, using the CCNR criterion. Ac- 
cording to Theorem 1, one can obtain an upper bound 
p = 0.9939, a = 0.232 for p{p) which is still entangled. 



FIG. 2: Figure of example 2 

When p = 0.9955, a = 0.236, using steepest descent 
method ,we get the concurrence is 0.101855, using our 
theorem ,the lower bound is 0.000487. Whenp = 0.9939, 
a = 0.232, using steepest descent method ,we get the con- 
currence is 0.101758, using our theorem ,the lower bound 
is 0.000019. 

Example 3. Consider the following states introduced 
inRef. 



FIG. 1: Figure of example 1 



Example 2. Pawel Horodecki introduced a 3 x 3 
bound entangled state in Ref. [28|, and the density ma- 



p{a) = -|* + )(*+| + -o+ + 
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-<7_, 



(15) 



where 2 < a < 5, a+ := ±(|0)|1)(0|(1| + |1)|2)(1|(2| 
|2)|0)(2|(0|), a := |(|1)|0)(1|(0| + |2)|1)(2|(1| 
|0)|2)(0|(2|). 
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FIG. 3: Figure of example 3 

In summary, we have provided an analytical formula 
for a lower bound of concurrence, by finding a connection 
with the currently most powerful detection criterion. The 
bound is very close to the actual values of concurrence for 
some special class of quantum states. Also, this bound is 
experimentally implemcntable and computationally very 
efficient, allowing to not only detect, but also to quantify 
entanglement in an experimental scenario. 

Acknowledgment. This work is supported by NSF of 
China(10901103), partially supported by a grant of sci- 
ence and technology commission of Shanghai Municipal- 
ity (STCSM, No. 09XD1402500). 



[1] R. Horodecki, P. Horodecki, M. Horodecki, K. Horodecki, 

Rev. Mod. Phys. 81, 865(2009). 
[2] O. Giihne, G. Toth, Phys. Rep. 474, 1(2009). 
[3] A. K. Ekert, Phys. Rev. Lett. 67, 661(1991). 
[4] C. H. Bennett, G. Brassard, C. Crepeau, R. Jozsa, 

A. Peres, and W. K. Wootters, Phys. Rev. Lett. 70, 

1895(1993). 

[5] R. Raussendorf and H. J. Briegel, Phys. Rev. Lett. 86, 
5188(2001). 

[6] W.K. Wootters, Phys. Rev. Lett. 80, 2245 (1998). 

[7] B.M.Terhal and K. G. H. Vollbrecht, Phys. Rev. Lett. 

85, 2625(2000). 
[8] K. G. H. Vollbrecht and R. F. Werner, Phys. Rev. A 64, 

062307(2001). 

[9] P.X. Chen, L.-M. Liang, C.-Z. Li, and M.-Q. Huang Phys. 
Lett. A 295 175(2002); E. Gerjuoy, Phys. Rev. A 67 
052308(2003). 

[10] A. Lozihski, A. Buchleitner, K. Zyczkowski, T. Wellens, 

Europhys. Lett. 62 168(2003). 
[11] K. Audenaert, F. Verstraete, and B. De Moor, Phys. Rev. 

A 64, 052304(2001). 
[12] F. Mintert, M. Kus, A. Buchleitner,Phys. Rev. Lett. 92, 

167902(2004). 

[13] F. Mintert, A.R.R.Carvalhoa, M. Kus, A. Buchleitner, 



Phys. Rep. 415, 207(2005). 
[14] K. Chen, S. Albeverio, S.-M. Fei, Phys. Rev. Lett. 95, 
040504(2005). 

[15] H.-P. Breuer, J. Phys. A: Math. Gen. 39, 11847(2006). 
[16] J. I. de Vicente, Phys. Rev. A 75, 052320(2007); ibid 77, 

039903(E) (2008). 
[17] S. P. Walborn, P. H. Souto Ribeiro, L. Davidovich, F. 

Mintert, A. Buchleitner, Nature 440, 1022(2006). 
[18] F. Mintert and A. Buchleitner, Phys. Rev. Lett. 98, 

140505(2007). 

[19] C. Schmid, N. Kiesel, W. Wieczorek, H. Weinfurter, 
F. Mintert, A. Buchleitner, Phys. Rev. Lett. 101, 
260505(2008). 

[20] Z. Ma, F. Zhang, D.-L. Deng, J-L Chen, Phys. Lett. A 

373, 1616(2009). 
[21] R. Augusiak, M. Lewenstein, Quant. Inf. Process 8, 

493(2009). 

[22] Z.-H.Ma,M.-L.Bao, Phys. Rev. A 82, 034305(2010). 
[23] I.Sargolzahi, S.Y.Mirafzali, M.Sarbishaei, Quantum Inf. 

Comput. 11, 0079(2011). 
[24] C.-J. Zhang, Y.-X. Gong, Y.-S. Zhang, G.-C. Guo, Phys. 

Rev. A 78, 042308(2008). 
[25] Y.-F. Huang, X.-L. Niu, Y.-X. Gong, J. Li, L. Peng, C.-J. 

Zhang, Y.-S. Zhang, G.-C. Guo,Phys. Rev. A 79, 052338 

(2009). 

[26] A. Peres, Phys. Rev. Lett. 77, 1413(1996). 

[27] M. Horodecki,P. Horodecki, R. Horodecki, Phys. Lett. A 
223,1(1996). 

[28] P. Horodecki, Phys. Lett. A 232, 333(1997). 

[29] K. Chen and L.-A. Wu, Quantum Inf. Comput. 3, 
193(2003);O. Rudolph, Physical Review A 67, 032312 
(2003); S. Albeverio, K. Chen, and S.M. Fei, Phys. Rev. 
A 68 062313 (2003); P. Wocjan and M. Horodecki, Open 
Syst. Inf. Dyn. 12, 331(2005); O. Rudolph, Quant. Inf. 
Process 4, 219(2005). 

[30] C.-J. Zhang, Y.-S. Zhang, S. Zhang, and G.-C. Guo, 
Phys. Rev. A 77, 060301(R)(2008); O. Gittsovich, 
O. Giihne, P. Hyllus, J. Eisert, Phys. Rev. A 78, 
052319(2008). 

[31] D. Brufi, J. Math. Phys. 43, 4237(2002);G. Toth, 
O. G-Uhne, Phys. Rev. Lett. 94, 060501(2005); D. 
Chruscinski and A. Kossakowski, Open Systems and Inf. 
Dynamics 14, 275(2007); D. Chruscinski and A. Kos- 
sakowski, J. Phys. A: Math. Theor. 41, 145301(2008);M. 
A. Jafarizadeh, N. Behzadi, Y. Akbari, Eur. Phys. J. D 
55, 197(2009). 

[32] O. Guhne, et al., Phys. Rev. A 66, 062305(2002); 

[33] M. Barbieri et al., Phys. Rev. Lett. 91, 227901(2003);M. 
Bourennane et al., Phys. Rev. Lett. 92, 087902(2004); H. 
HaiTner et al., Nature 438, 643 (2005). 

[34] K.Chen and L.A. Wu, Phys. Rev. A 69, 022312 (2004). 

[35] R.A. Horn and C.R. Johnson, Topics in Matrix Analysis 
(Cambridge University Press, New York, 1991). 

[36] P. Rungta, C. M. Caves.Phys. Rev. A 67, 012307(2003). 

[37] C.H. Bennett et al., Phys. Rev. Lett. 82,5385 (1999). 

[38] S. Ryu, W. Cai, and A. Caro, Phys. Rev. A 77, 
052312(2008). 

[39] M. Horodecki, P. Horodecki, R. Horodecki, Quantum In- 
formation An Introduction to Basic Theoretical Concepts 
and Experiments Springer Tracts in Modern Physics 
(Springer, Berlin), (2001). 



